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sf ' ABSTRACT u)b^d ^ 

Ac reques t ed , we have examined the Venus Ballooi\ concept in 
order to further develop the ideas and concepts behind it, and to 
creatively apply them to the design of the major Venus Balloon 
components. This report presents our models of the vertical path 
taken by the Venus Balloon and the entry into Venusian atmosphere. 
It also details our designs of the balloon, gondola, heat exchanger, 
power generator, and entry module. A vehicle is designed for a 
ballistic entry into the Venusian atmosphere, and an atmosphe ric 
model is created. The model is then used to set condi tions^The^ shape 
and material of the vehicle are optimized, and the dimensions of the 
vehicle are then determined. Equipment is chosen and detailed that 
will be needed to collect and transmit information and control the 
mission. A gondola is designed that will enable this sensitive 
electronic equipment to survive in an atmosphere of very high 
temperature and pressure. This shape and the material of the shell 
are optimized, and the size is minimized. Insulation and supporting 
structures are designed to protect the payload equipment and to 
minimize mass. A method of cooling the gondola at upper altitudes 
was established. Power needs of the gondola equipment are 
determined. Power generation options are discussed and two 
seperate thermo-electric generation models are outlined. 
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VENUS CLOUD BOBBER MISSION PROPOSAL 


SYSTEM DESCRIPTION 


An entry vehicle and survey package are proposed for a 
scientific expedition to Venus. The total mass of the vehicle and 
package is 157 kilograms (345 pounds) and has a total volume of 
0.179 cubic meters (11040 cubic inches). The entry vehicle will be 
attached to a spacecraft bus that could possibly carry multiple probes 
of this same configuration. The entry vehicle is a trimmed blunted 
cone of half-angle 45 degrees with a base radius of 19 inches and a 
bluntness ratio of 1/2. An AVCOAT 5026 heat shield of thickness 2.81 
inches protects a titanium aeroshell of 0.164 inch thickness. Inside the 
entry vehicle is a survey package with a total mass of 89 kilograms 
(196 pounds) which consists of a dual balloon system, a science 
gondola, and a heat exchanger. The cigar shaped polyethylene 
balloons carry hydrogen as the primary gas and R30 as the phase 
change material for a cylindrical aluminum heat exchanger. The 
balloons have a mean radius of 2 meters (7j> inches), thickness 4 milK 
and combined height of 22 meters (866 iSctfes). The heat exchanger is 
0.3117 meters (12.27 inches) long and the diameter is 0.2286 meters (9 
inches). It serves as altitude control for this survey package, as an j/r 

unmodified dual balloon system would reach an equilibrium height.-'* 

The maximum range of the balloon system with the heat exchanger 
is 64 kilometers and is attained three hours after deployment from 
the entry vehicle. The gondola is a single titanium sphere of radius 8 
inches and thickness 90 mi]$Tlt has fibrous aluminum-silica 
insulation on the inside of thickness 1.41 inches. Power is generated 
by a smaller scale heat exchanger, and is used to run the camera, 
sensors, and communications equipment onboard the gondola. 


Stucco's seen.******* 
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MISSION OBJECTIVES 


The primary goal of the mission is long term observation and 
measurement of the low altitude properties of the Venusian 
atmosphere and terrain. Other possible science returns include the 
study of wind variation by time of day and season, evaluation of 
gondola decay due to acidic exposure during ascent and descent 
through the cloud layer, and a precedent setting d esign that could be 
used for surveys of all planetary atmosphere^ Adaptations of the 
basic gondola desigrvwould allow surface sampling missions and 
passive horizontal flight control. The mission is desirable because of 
its lightweight design, which reduces payload costs. The mission also 
offers the prospect of a repeatable design, where several of these 
survey packages could be produced, reducing the high cost of 
producing a unique one shot package. 
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The simplest method of entry was chosen for study. The 
ballistic entry is a trim, minimal heat trajectory characterized by a 
constant angle of entry, four basic equations of motion, and good 
closed form solvability (Weisel). It is required that the vehicle have 
negligible or zero lift, or skipping may occur. A skip entry is 
characterized by an entry vehicle that pulls out of its descent and 
exits the atmosphere at the same angle it entered (Regan). The craft 
eventually descends again, with reduced speed. The skip entry would 
be a considerable alternative if peak heating rates for ballistic entry 
were found to be too great a burden for a lightweight vehicle to 
handle (Ely). Moderate lift to drag ratios produce partial skipping, in 
which the craft momentarily ascends without leaving the 
atmosphere before falling to the planet (Regan). 
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Deployment Sequence 

In all likelihood , the decelerate 
spin rate. There are some advante ^ ^ 

unchecked. A slow, steady spirt wiU caiK^ out)u^xpected lift forces -v 
and help distribute thermal energy more evenly tDuncan). Currently/ 
no provisions exist for inertial navigation or propulsion systems. 

However, a simple yo-yo d^pin mechanism could be used to set the 
spin rate before entry (Weisel). Such amechanism was used with the 
Pioneer probes and proved to be ^ycheapmeans of spin axis control - 
(Fimmel). Even though the vehicle is supposed to b^trimmed, any 
uncompensated lift could increase entry time, resulting in increased 
heating and less accurate tracking of the vehicle as it descends. To 
prevent tumbling of the craft during entry, the center of pressure 
should be located extremely close to the center of mass (Griffin). The 
desired orientation of the entry vehicle is shown in Figure 1; note 
that the only major requirements for attitude control are 
maintaining colinearity between the craft's axis of symmetry and 
keeping the craft in the nose-first position instead of the base-first , 

position. A balance will have to be made between the desire to push "7 Oic&D < 
the mass center towards the nose of the vehicle, and the need to keep [ 
thermally sensitive components away from the hottest portions of 
the craft. 
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Figure 1: Forces on a ballistic entry 
vehicle (Weisel, p.220) 

Equations of Motion 



m 


dX 

— = V cos (angle) 
dt 


dV 

mV — = -mg cos (angle) 
dt 


dH 

— = V sin (angle) 
dt 


dV 

m — = -D - mg sin (angle) 
dt 
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These equations of motion can be manipulated to solve for the 
velocity (V} z _5iidJhe_aerodynamic deceleration (a). For each height 
listed in<^e^readshee^hereare five unknowns: ^ 

Kd/ an^ w e rs e~dmsit}rpafameter ■ 

Ho, a fi^ious reference height jrr mx ^er *esv**s^ 

Bo, an energy parameter ^s^rr/o^aa 

a, the aerodynamic deceleration, and 
V, the craft's current velocity. 

^ ‘"’"""A complete list of terms used in this report is available in 
appendix A. Five equations are available to solve for these five 
unknowns, and ajn^re-detaiied treatment of determining the height 
profile is given irf^\ppendix ^ _ 
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The deployment sequence begins when the craft speed has 
decreased to 35 meters per second. An accelerometer and pressure 
altimeter will be used to determine when this speed has been 
reached. At that time the onboard computer will send a command to 
release the chute attached to the aft cover. This will further decrease 
the craft speed. When the craft speed has been reduced to 10 meters 
per second, the computer will activate the explosive bolts and inflate 
the primary balloon. Sensors in the canister valves will monitor the 
flow of gas into the primary balloon. When the computer receives 
readings of negligible gas flow, the computer will send its final 
command to sever the survey package from the aeroshell. The 
sequence is shown in Figures 2a-2d. The entry vehicle has served its 
purpose and the operational phase of the survey package can now 
begin. The final speed at which the survey package begins the 
operational phase is about 9 meters per second. 





A Model Atmosphere 


/?r 


Data on the physical composition of the Venusian atmosphere 
was obtained for heights up to 120 km by Pioneer probes (Fimmel, 
Vargaftik). Above 120 km, a mathematical model was used (Hunten). 
This model interfaces smoothly with Pioneer data and serves as a 
reasonable approximation. The variations of temperature, density, 
and gravity with height for this model are shown in Figures 3a-3c. 
Instead of taking a conventional simplified treatment of the 
atmosphere for entry calculations, the data was used to create an 
accurate model of vehicle descent. The traditional entry model uses 
an isothermal exponential atmosphere of invariant composition 
(Weisel). While this makes the determination of the trajectory and 
heating profile easier, it sacrifices a significant amount of accuracy. 




Figure 2a: Initial deployment 
of parachute 


Figure 2b: Velocity reduced 

from 35 m/s to 10 m/s 



Figure 2c: Explosive bolts activated 
and aft cover released 


Figure 2d: Balloon inflated, 

connections severed 
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Figure 3c: Variation of Gravity with Height 



The density graph is similar to an exponential atmosphere. The 
gravity graph has nearly a linear look, since the distance covered is 
much smaller than the planetary radius. The temperature graph is 
disjointed somewhat at the spots where the mathematical model and 
the probe data are supposed to mesh, but it does not have a 
significant effect on the deceleration and heatng profiles. 
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The model used to calculate deceleration and heating rates has 
values for pressure, density, temperature, gravity, viscosity, and 
molecular weight of the Venusian atmosphere in discrete height 
steps that traverse the lower 200 km of the atmosphere. In reality, 
the atmosphere extends upwards of 500 km; significant deceleration ^ ^ 
and heating do not occur above 200 km for this vehicle The ** 

atmospheric properties allow calculation of a deceleration profile 
with respect to height; the velocity and time are found by integration 
assuming a step deceleration function. The velocity profile enables 
calculation of the convective heat rate from the atmosphere to the 
craft (Dueber). The time profile allows the use of an energy method 
in determining the necessary shield and aeroshell thicknesses. 

Sensors on the skin of the vehicle could provide valuable data to help 
confirm current models of the Venusian atmosphere, as well as 
increase understanding of the entry heating problem. 

Shape Selection 
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Four basic craft shapes were considered and are shown in 
Figures 4a-4d. The shape of the craft is important, since it directly 
influences the deceleration and heating profiles. A trend emerged 
from comparisons of the different shapes at a set volume: truncated 
and blunted cones produce significantly lower heating rates and 
lower amounts of total heat that the craft must deal with. Between 
these two shapes, four elements determined the final shape selection./ Zeas ^ 
The blunted cone had smoother heating and deceleration profiles due/ ^ ^ * 
to a higher drag and a lower area. The truncated cone had reduce d / * ^ 
entry times and could endure greater entry angles due to a smaller J 
drag a nd an increased area. The fin al decision came down to the 
preference of red uced peak heating rates over reduced entry ^ *s~s\ 

addition to t he minimization of the necessary aeroshell and shiPlH f ^ ***** 
thicknesses,. The blunted cone became the chosen shape of the vehicle, [ roK'c r&rz- J 
but now materials and dimensions needed to be optimized for the '*’**’• / 

package that must be protected inside the vehicle. See Appendix C 
for the numerical analysis of the four shell shapes. " 
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Figure 4a-4d: The four shapes considered 
for the entry vehicle 

At the time of design, the dimensions of many internal 
components were tenuous at best. Other criteria were usedjhowever, 
to select optimal cone half angles and bluntness ratios. A high draj 7 
coefficient is desirable to ensure that the vehicle decelerates to an 
acceptable speed in a limited amount of time. A large surface area 
can also aid deceleration, with a penalty of increased heat transfer. 
The entry constraints were assumed to be a maximum decelefc 
of 500 g's and an entry time close to 5 minutes. Four size iterations 
produced the optimal cone angle of 45 degrees and the bluntness 
ratio of one half. This is the same basic_sh a Pe as the Pioneer small 
probe, with only a difference in sec 






Material Selection 


Two important aspects of the vehicle needed to be determined. 

The heat shield and aeroshell of the vehicle will be subjected to the ^ ^ / 
most adverse conditions, and must be designed properly to in sure— — • 
protection of the survey package. The heatshTeidancTshell must be 
able to handle high dynamic and thermal loadings. The heat shield 
ideally has the properties of low density, low thermal conductivity, 
high heat of vaporization, and high emissivity (Dueber). Whether 
AVCOAT was charred or virgin material, it ranked higher than other 
candidates in almost all of these properties. The only directly 
determinable input is the convective heating from the atmosphere to 
the craft (Regan). This input can be converted to a total energy 
transfer to the vehicle. An approximation was made that the 
required mass of ablative heat shield material equals the total energy 
transferred, divided by the heat of vaporization for the shield. After 
finding that the mass required was probably too small, a more 
conservative approach was taken. The initial and fuia| Jdneticx^ 
energies of the vehicle were calculated, and und efTPU advice the 
heat shield mass was taken to be approximately 20% of the total 
craft mass. Under these conditions, the craft will empty 15% of the 
kinetic energy through mass loss. The remaining 85% will have to be 
radiated away or conducted through the aeroshell. Appendix C has „ _ 
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the numerical results of shell and shield material studies. 


Dimensioning the Entry Vehicle 






The first obstacle in finding the optimal dimensions of an entry 
vehicle is the establishment of a complete equipment list. Once this 
list has been made to the best possible accuracy, the total mass and 
volume are used to set up a 500g entry. Ihg 

spreadsheet returns values of the maximum convective heat flux, the 
convective energy seen by the craft, and the angle at which the entry 
occurs. These values are used in heat transfer and stress calculations 
to determine if the shield and aeroshell thicknesses are accurate. The 
process is qfcfte iterative, and time consuming unless computer aids 
are used. A detailed explanation of thgjieai_ri3m>fer and stress 
calculations are found in Append icef^U and E^^he current equipment 


list is shown on Table 1 and the spacecraft dimensions are shown in 
Figure 5 ' / — ■ — 
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Totals Allocated 1 57 kg (345 Ibm) 0. 1 79 m3 n 1 040 |n3t 

Craft design based on 150 kg (330 Ibm) 0.180 m3 (11100in3) 
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half cone angle 



tshield, shield thickness 


SPACECRAFT DIMENSIONS 

Base Radius, R: 0.482 m (19.0 in) 

Nose Radius, Rb: 0.241 m (9.49 in) 

Half Cone Angle: 45 degrees 

Shell thickness, tshell: 0.416 cm (0.164 in) 

Shield thickness, tshield: 7.14 cm (2.81 in) 
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Figure 5: Exterior dimensions of the vehicle 
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The aeroshell is the next line of defens e/^ince temperature 
distribution through the thickness is unknown, the criteria that drive 
shell selection are basic: Low coefficients/of thermal expansion, good 
buckling resistance, low oxidati©»=^ates> low thermal conductivity, 
and lightweight materiaT^Jaworski)^ worst case heat transfer 
scenario was developed^Trsirtg-tb^peak convective transfer rate as a 
steady state input to the exterior of the shell. The interior 
temperature was set at the desired level, while the exterior was //. 
treated as a hot spot at an arbitrary high temperature.^t the time , "2 
this seemed^dre me, but it turns out that thp Hynamir lnark J 

the shelTto be even thicker than calculated here . Only two of four 
material candidates had desirable properties; between titanium and 
stainless steel, titanium handles large pressure vessel stresses with 
lightweight, thin sections. Titanium was chosen to be the aeroshell 
material, and was designed to sustain a deceleration induced hoop 
stress of 0.8 yield and handle thermal hoop stress of 025yueld 
(FortescueJTSince peak deceleration and peak heating do not occur 


^0 


at the same time, it is reasonable to assume that no combination of 


the two stresses causes failure. 


£Zrr G&V Cctyeter 


//z SC> ^7^ 


It is regrettable that the depletion profile of the heat shield OrfjCCA*# yy/l/C* Jy 

could not be established. It leaves some doubts about the safety of the 
survey package should the thermal conditions prove to be worse 
than anticipated. This drove the decision to adhere a thermal blanket 
as a f inal line of defense , (see Figure 6) This is the same thermal 
blanket that was used on the Magellan orbiter. From the figure, it is 
noted that the thickness of the blanket is greater than the aeroshell 
thickness. It is a small price to pay for added thermal protection. 

Originally, the interior of the vehicle had been partitioned by shelves; 
a(QPL^ehest these have been removed and it is assumed that there 
vdH be attachment points on the inside of the shell for the survey 
f \ package and the rest of the deployment equipment. 
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Figure 6: Suggested thermal blanket 

(Miller, p.8) — y a^ra> c** **-- 

Assumptions: A General Disclaimer 

Entry calculations and the design of an entry vehicle are 
invariably a very complicated problem. It is important to realize the 
many simplifications and assumptions made during the design 
process, as the accuracy and reliability of the design will be directly 
dependent upon them. Let's begin the discussion by examining the 
accuracy of the atmospheric model. The most important portions of 
the model are dependent upon the accuracy of Pioneer sounder probe 
data. It would be prudent to expect cyclical variations in the physical 
characteristics of the atmosphere (Hunten). The magnitude of these 
variations, while unknown, could have a sizable influence on the 
heating and deceleration profiles in the middle altitudes. Gravity has 



been approximated using the inverse square law for a surface value 
of the Venus gravitational acceleration, and the angle of entry has 
been assumed to be constant for the entire entry. The spacecraft mass 
and cross sectional area have been assumed to be constant during 
entry. In reality, neither the angle, mass, nor area will remain 
constant. The drag coefficient, which should vary over the entire 
trajectory, has instead been taken as an average value over the 
entire flight. The time elapsed is based on integrating stepped i 

deceleration values between heights. Convective heating is based on ^ , 

an average skin coefficient for the entire vehicle, although this ^ ^ - 

parameter is allowed to vary with respect to height. The viscosity of ^ 
the atmosphere is linearly interpolated at heights of 0-50 km, and / 
assumed to be ’ J ’ ' x 
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from 51-200 ki* (Vargaftik). . . . 




Sensitivity Overview 


The design of the craft suffered with respect to thermal 
modelling. The true energy balance for the craft should look like 
Figure 7. Although factors relating to the mass flux can be calculated, 
the mass flux itself is a function of vehicle wall temperature (Kays). 
The problem is essentially a fourth order nonlinear nonhomogeneous 
transient partial differential equation. If it sounds diffirnlt try — 
solving it and find out what the designer faced . The energy lost to 
radiative emission should reduce the conductive energy transferred 
through the wall. The design has an ample amount of shielding for 
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design conditions where ablative dissipation accounts for up to 15% • 

of the total dissipation requirement. Beyond that, hopefully. thp^ -^ ] 

aeroshell and thermal blanket can keep the interior temperature at a ■ 

safe level. The use of a temperature and composition variant 
atmosphere increased the accuracy of the deceleration and heating 
profiles. Only further testing of the Venusian atmosphere and direct 
recording of flight data can increase the accuracy of that part of the 
design. The use of pressure vessel analogies to quantify the stress 
behavior of the aeroshell may come into question: worst case 
scenarios involve limited yielding of the titanium aeroshell with an 
acceptably low risk of failure. The shearing of the shell has not been 
directly considered, as the dynamic loading is a distributed one and 
does not solve easily for this shape. T he location of components, 
which was not readily established, could greatly affprt thp dvnamir 
stability of the craft. — 
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Ablation Energy 


q - Convective Flux Ablation Energy = mfhvap 

A*T A 4 = Radiative Emission Conduction = -k*(dT/dx) 

where mf is the mass flow of the heat shield 
hvap is the heat capacity of the wall 
k is the thermal conductivity of the wall 
dT/dx is the thermal gradient across the wall 
A is a constant based on emissivity 

Rate Balance: -k*(dT/dx) = q" - A*T A 4 - mf*hvap 

where q'' and mf are unknown functions of T 


Figure 7: The energy fl ux ba 
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Expansions of Work 
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While much has been done, many aspects of vehicle design have 
not been covered completely. The thermal balance of the spacecraft 
needs refinement. Layout and interfacing of equipment needs to be 
developed. Calculations of dynamic and thermal loadings have been 
greatly simplified and warrant closer inspection. The design of the 
entry vehicle can really be broken down into four sections: aeroshell 
design, heat shield design, selection of subsystems, and internal 
layout. Each one of these subjects is interrelated and would be best 
addressed by a one year developmental effort. It is hoped that 
further investigative work will be done by future students. 






Knowledge Gained 
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In learning how to design an entry vehicle, the designer gained 
some proficiency in certain topics. The creation of an accurate 
atmospheric model was undertaken successfully, as was an 
alternative material selection process in the absence of some useful 
d^a. The theory behind heat shield selection and ablation physic^ 
was understood, if not successfully applied. Conversations witltfPlTy 
revealed the tendency for designers to stick with existing shapesln 
vehicular design. A significant improvement was made in the 
acquisition and use of resources for research. The truth is that more 
time and experience will produce a better design, but the effort 
described here is by no means a token one.^ ^ 

Transitional 


The entry vehicle requires little or no power during descent and 
deployment phases. During descent, atmospheric data will be sent to 
instruments onboard the gondola for storage. With the deployment 
completed, the next design obstacle is a heat exchanger that will 
regulate the survey package’s altitude. The development of a 
gondola package is considered next. 

7T Gondola ^ qt ") 


The gondola for the V enus balloon is a complex problem, 
because of the necessity of keeping temperatures inside the gondola 
below 20° Celsius while the surface temperature of Venus is 460° C. 
This is necessary because the electronic equipment within the 
gondola cannot be subjected to cyclic high temperatures without 
sustaining damage. Also, the gondola must be strong enough to 
withstand a pressure of up to 92 atm. 

The first steps in designing the gondola were to estimate the 
size and to decide upon a shape. A list of equipment for the 
gondola is shown in Table 2. (^ffie jlhen obtained volume and mass 
data for the chosen equipm^ntSome of the data^dbtained was 
estimated or abridged ir/order to obtain what^felt were more 
accurate numbers. Based on these dimension values, a total payload 
area was found to be 726.6 in? . At this point, ^e^ecided that due to 
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GONDOLA MASS ESTIMATE 


COMPONENT 

MASS 

Power Distributor 

3.5 kg 

Charge Regulator 

3.0 kg 

Battery 

3.3 kg 

Antenna 

0.7 kg 

Transmitter 

3.5 kg 

Camera 

2.0 kg 

Command and Data Storage 

2.4 kg 

Shell 

5.234 kg 

Insulation 

0.6902 kg 

Foam 

.3195 kg 

PCM and Heat Pipes 

1.5 kg 

Miscellaneous 

0.5 kg 

TOTAL MASS = 

26.64 kg 
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Tgfcle 2. Gondola Equipment Mass Estimates 


it’s simplicity, aerodynamic features, and capabilities as a pressure 
vessel, a sphere would be the best shape to make the gondola/with "> 

this in mind/^vCXv orked backwards fromc^up 1 payload volume 

estimate to find that the sphere should have a minimum fadius of 5.6 

in. To ma ke ro om for insulation, support structure, and other added f 

v olumesxdypselection for the sphere size was a radius of 8.0 in.- - r ^ 

The next step in our design process was to select the material 
and thickness of the gondola. Due to it’s excellent strength to weight n 

ratio a titanium alloy (6% Al, 4% V) was chosen.'Titanhnn resists^ 
hr struts and thin plates be tter than other common 
JpattTrials, such as stainless steel and aluminum. Titanium also has a 
yield point for axial stress. Based on these material properties, 
the atmospheric pressure of Venus, and the thin-walled pressure 

vessel formulas, a thickness of 0.04508 i n. was found in hg criti cal 

value. In order to protect against atmospheric pressure variations, ^ 

possible internal point loads, and thermal deformations, a safety 
factor of 2 was used to bring the thickness of the sphere to 0.09016 in. 

Based on this number and the density of the material, a mass value of 
5.234 kg was obtained for the sphere. 
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At this point in our design, there were two 
were being considered. One type consists 
with a vacuum in between as an insulat 

insulation would be very effective at preventing both conduction and 
convection, and also would be very lightweight. One problem with it 
is the fact that the two spheres would have to be attached in order to 
provide structural stability and support. ^Using these supports- — " 
undermines the conduction resistant properttes-of^e\^a^um, as 
well as applies unwanted point loads to each sphere. It also 
maximizes the pressure difference acting on the sphere because one 
side of the sphere is at approximately 92 atm and the other side (the 
vacuum) is at 0.0 atm. In addition, having two shells would increase 
the weight of the gondola considerably. 

The second structure considered was^Tto use one spherical 
shell, and line it with enough insulation so that the heat transfer to 
the equipment would be minimal. With this type of structure there 
would only be one shell, so there would be a weight decrease of the 
titanium used. In addition, this type of gondola design makes 
supporting structures much simpler and more efficient. The main 
drawback is that there is a weight gain due to the insulation used. 
Based on the advantages and disadvantages of these two designs, 
we chose to pursue the single-sphere model. A 

Insulation 

The next step in the design process was to determine thetype 
and thickness of the insulation material. The driving factor was the 
weight vs. the insulating ability, but a high maximum service 
temperature was also required. Based on these factors, we chose to 
use a n alumina-silica fiber as the insulating materia l^. Us i ng fhe 
material properties of this fiber, the temperature differencefrBiri the 
Venusian atmosphere, and the desired operating temperature inside 
the gondola of 15° Celsius, it was calculated that the thickness of the 
insulation would have to be 1.41 im Using this number and the 

density of 48 kg/m 3 , the masspHhe insulator is 0.6902 kg (Incropera 
and Dewitt). This proved tp45e much less than the mass of the 
additional shell which wofild have been needed to provide vacuum 
insulation, so the chokfe to discount the vacuum insulation seems 
reasonable. 
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Support 


The next part of the problem was to design some sort of 
supporting structure to hold the equipment in place. In addition to 
support, other needs were to absorb shock and dampen vibratimrg. 
/^We -chose to use a molded polystyrene foam material for this.<WfD 
decided to envelope the equipment and phase changing material 
completely in foam so that the entire otherwise unoccupied area 
inside the shell would be polystyrene. This foam was of suitable — ' 
* stiffness and is operable in somewhat higher temperatures - — 
(Klempner and Frisch, p. 175). This serves the needs stated above as 
well as adds another insulating layer. After subtracting the volume 
of the equipment, phase change material, insulation, and all other 
components from the available space inside the gondola^v^found 

that there were still .019966 m3 of volume left over. This volume of 
the chosen foam has a mass of .3195 kg. 


Cooling 


Despite having sufficient insulation from outside heat, it is 
important not to forget about the considerable amount of heat that 
the electronic equipment itself can generate. Because of this, we also 
needjomdude some type of cooling device, (^©decided to use^fv j 
'5ctadecane)as the phase changing material, along with a pair of 

rediode reflux heat pipes. The phase change material (pcm) 
will be contained in two lightweight heat conducting plastic 
containers and placed around parts of the equipment. The pcm will 
absorb heat from the equipment. When the temperature of the 
equipment reaches the melting point of the pcm, it will melt and by 
doing so absorb a greater amount of the heat energy. The heat pipes 
will be used to cool off the phase change material (pcm) when the 
balloon is in the upper atmosphere. They will conduct heat from the 
pcm through the foam and insulation and onto the shell, where it will 
convect into the atmosphere. This heat transfer will continue until ft 
the pcm has reached a much lower temperature and has re-frozen 
The mass of the pcm and heat pipes is estimated at 1.5 kg. [j 
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Extra Mass (Equipment) 


In addition to all of the previously detailed equipment, there 
are also several other sources of mass and volume to consider. 

There will be weight for electrical wires to connect each piece of 
equipment. There will also be a need for a small length of fiber optics 
cabling to allow the camera to take pictures outside of the gondola. 
These additional masses are estimated to be 1.0 kg. 

Power Generation 



} 


Due to the large power requirement of the Venus balloon and 
the small allowable mass of the gondola, power must somehow be 
generated on Venus during the mission. The initial step taken was to 
determine how much energy would be needed to operate the balloon 
during one cycle. The equipment is already known and also the 
power needed to operate each piece of equipment. This data is 
shown in Table 3. As shown in this table, the estimated power 
requirements of each piece of equipment are multiplied by the 
estimated operating times. These estimations are based on the 
length of one cycle of the balloon being 24 hours long. This gives the 
energy requirement per cycle for each piece of equipment. These 
numbers were summed to find the total energy requirement of the 
gondola per cycle, which is 908.7 kj. —z? 


GONDOLA EQUIPMENT 
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POWER 


OPERATION TIME PER CYCLE ENERGY USE PER CYCLE 


0.5 W 

4.0 W 
N/A 
N/A 
40 W 

5.0 W 
2.4 W 


86400 sec 

43200 J 

3600 sec 

144000 , 

7200 sec 

288000 , 

300 sec 

1500 J 

86400 sec 

432000 , 

TOTAL ENERGY PER CYCLE = 

908.7 kJ 


Table 3. Gondola Equipment Power Consumption 





The next step m^oupdesign was to find a good method of 
power generation, (^elooked into using solar power, wind power, 
and thermal power. Solar power was the first option looked into. 

Although it is often used in space applications, that is typically on 
orbitting sattelites. It was discarded fairly quickly for a number of 
reasons. For instance, solar arrays, despite continuing 

improvements, have low output. They are typically 12-15% efficient / pat& 
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(Boer). Also, the cloudy atmosphere, combined with the facts that it 
will be night for roughly half of each cycle and that there is a loss of 
solar radiation flux for a tilted solar panel, means that the solar 
intensity, and therefore the power output will be very minimal and 
sporadic (Rapp, p. 37). In addition, high temperatures cause damage 
to solar arrays and loss of efficiency. 

Turbine power was thought to be an option for use due to the 
fact that the gondola would have constant vertical motion. The 
- vertical air flow caused bv this motion would then be able to drive 
vertical wind turbines.(J©did not look closely at the power 
capabilities of this option.<®>decided that we did not want to have 
a moving mechanical part like a turbine due to the vibrations and 
cyclic stresses that would be present, affecting both the turbine and 
the gondola. It would add difficulty to the design of other 
components, such as the gondola shell and the equipment. In 
addition the vertical velocity of the wind could often be as low as 1 
m/s, which would not seem to offer a great supply of energy. 

Finally, using wind turbines would add extra weight by requiring the 
use of a generator. 

Therefore, $ve^:oncluded that $v£^hould try to harness some of 
the vast thermal energy found on the planet Venus. The methodctvg^ 
decided to use was a thermocouple. The thermocouple device shown j 

i i. , ,, ... . - j 


can be used to create a hot and cold node to which a thermocouple 
can be attached. The thermocouple itself has not been designed. 
although an efficiency value of 6% was assumed. Although an 
efficiency value of 10-20% is possible in a tljgrmocouple, the chosen ^ c^a / V' aa£>7~ 
value to use in calculations is taken to be(6%^due to the harsh 
operating conditions in the atmosphere of Venus. Using this 
efficiency , it can be calculated that 1 6333 k/ of heat energy would 
need to be converted into electrical energy in order to sufficiently 
power the balloon. A summary of the w^rk accomplished in 
designing this method follows. 






y 

a ^ 








^ _ ~frrcr a'-uzA 


wt*- ■ 

The model showjaTs a device that can be used to store this heat 
energy while it is-bdng converted Into electrical energy. Exact 
numbers have not yetbeen workedjout, but this analysis shows the 
scope of the problemT~Water has beenchosen as a working fluid, 
although there may be much better options. _ 

In stage one, as shown x the fluid is contained in section A. As 
the water is heatedupT lT'evaporates and raises into an adiabatic 


Lo 


section, B. Once this happens, the valve is closed and the balloon 
rises to the upper atmosphere. It has been estimated that an average 
temperature of 350 Celsius will be maintained in section B. That 


would give the water vapor a specific energy of 2877.4 kj/kg (Moran 
and Shapiro p. 700). A thermocouple is connected between section B 
an d an arbitrary cold poin t, which can be basically anywhere in the 
^upper atmosphere, perhaps on the gondola. At this point the balloon 
may oscillate in the upper atmosphere. According to this model, the 
majority of heat lost would be through the thermocouple and 
approximately 6% of that would be converted to electrical energy. If 
the heat energy loss were to occur over time until the water vapor 
became a saturated liquid at 100 C, where latent energy of 
vaporization could be taken advantage of, it would prove to be 
beneficial. At this point the specific energy of the water would be 
418.94 kj/kg (Moran and Shapiro p. 700). This is a difference of 
2458.46 kj/kg between the two states. Based on this number and the 
total heat energy needed, 6.16 kg of water would be needed to insure 
proper power generation. This is probably not a very reliable 
number, but th^re are many things which can be done to improve 
upon it. t - 

Many assumptions have been made due to lack of availability 
of accurate heat transfer modelling through the adiabatic section of 
the thermocouple device. Therefore there are many ways to improve 
this model. The working fluid and pressure can be changed to take 
greater advantage of specific internal heat energy. Also, the process 
can be reversible to some extent when the balloon descends. The 
adiabatic section will keep the water cooler than the surrounding 
atmosphere until the valve is open, so the thermocouple will have a 
temperature difference. Possibly the efficiency of the thermocouple 
can be approved. It is also possible to link the thermocouple directly 
to the heat exchanger of the secondary balloon. This would save 
weight for the overall design. It is too soon to predict how that 
would affect the power output of the thermocouple. 



• a In ^? nclu ^ lon ' having previously looked into solar arrays and 
wind turbines, the thermocouple seems to be the most promising 
choice for power generation and it warrants further study. The heat 
exchanger that will control the altitude of the balloon is now 
discussed. , „ 
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INTRODUCTION TO THE VERTICAL MODEL 


Objectives 


The purpose of this section of the project is to investigate the 
^P° rn l ance proposed Venus Balloon for exploring the deep 
°f the planet. Knowledge of the vertical motions of the 
balloon system is important for three reasons. First, to effect a 
desired ascent profile, we need to understand how tfie balloon 
^ s P°P ds to its lift force, drag force and changes in its operational 
environment. Second, to initiate ascent or to vary the rale, we must 
be able to predict the effects of gas vagjng. y 

to provide adequate cooling for the onboard instruments 
, o f F" lcient communication time, the altitude versus time plots 
must be carefully and accurately determined. * 


Background and Challenges 


P 16 temperature Venus atmosphere of mostly carbon dioxide 
poses a great challenge to reconnaissance work on the planet. 
Intensive and prolonged exploration of the planet surface and 
atmosphere have been impeded by the unfriendly atmosphere of 
j' The only balloons ever flown in the planet were the fixed 
altitude Soviet Vega Balloons which lasted for about t£o days 

The opacity of the thick cloud layer; 48 km to 60 km above ground 
level (a.g.l.J, and the haze layer; 3u km to 48 km a.g.l., precludes anv 
clear imaging from a higher altitude where the temperature and Y 
pressure is more moderate. The dense atmosphere also precludes 
the transmission of signals to Earth based stations (see Fig. 8a ) In 
order to obtain clear images of the terrain, the balloon has to 
descend to altitudes where temperatures of 460 Celsius and 
pressures of 92 bars prevail. ^ 1 

To communicate, the balloon has to ascend above the cloud layer. 



VENUS ENVIRONMENT 
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Sources of Information 


The first step towards achieving'lSurjobjectives was to review 
the rather extensive literature in one o HKree problem areas that 
deal with balloon ascent. The first area is the Venus atmospheric 
conditions. Sources consulted include: 


o Hunter et al., 1983 
o Fimmel, 1983 
o Lide, 1992 
o Nock and Jones, 1993 
o JPL correspondence: 

JPL scientists. Jay Wu and Jack Jones 


The second area deals with balloon ascent, weather balloon 
trajectory modelling, space exploration vehicles, float motion and 
aerodynamics of free pressure balloons: 




o Soviet Vep Balloon Report, 1982 

o Dwyer, 1985 
o Ward and Kincaid, 1985 
o Boyce, 1986 
o Zdk, 1988 
o Greenberg, 1988 


UW-Faculty consulted with: 

Professor Fieriesen, Fluid Mechanics, M.E. Dept. 
Professor R. Reitz, Fluid Mechanics, M.E. DepT. 
Professor G. Meyers, Thermodynamics, M.E. Dept. 
Professor P. Cheng, Thin Shells, E.M. Dept. 

ohnson. Mechanics, E.M. Dept. 

. Kuelbs, System of Eqns. Math Dept. 
Frosteric, Diff. Equations Math Dept. 


-i. A V/A 

Professor 

Professor 

Professor 


o JPL correspondence: 

JPL scientist. Jay Wu (Trajectory) 


The third area, heat transfer (internal and external) 
encompasses some of the principal concerns that compelled this 
stud y The altitude control of the balloon is directly related to the 
overall heat exchange rate of the system. J 
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C , Of the three problem areas, this is by far, the most complex, 
jt hout p erforming actual experiments, most heat exchange 
_,_jfH 3 ieter^cannot be determined to any level of -fair-accuracy. The 
closest thing to performing an actual experiment on expensive 
prototypes is computer simulation. Fluent, a commercial fluid 
dynamics package was used to investigate the heat exchange rate 
and the bulk temperature of the balloon as a function of time and 
altitude. Sources on heat exchange consulted include: 


o Seigel et al., 1981 
o Vargaftik, 1983 
o Dwyer, 1985 
o Blamont et al., 1985 
o Incropera et al., 1990 
o Fluent Manuals., 1993 

o UW-Faculty consulted with: 

Professor Mitchell, Heat Transfer, M.E. Dept. 
Professor Meyers, Thermodynamics. M.E. Dept. 
Professor R. Reitz, Fluent, M.E. Dept. 

Ms. C. Maul (T.A.), Fluid Mechanics Cnem.E.Dept. 



Scope of Investigation 


m 


This project involved seeking information on /designing for the 
following: 

o Balloon design 

Selection of working fluids for the balloon system 
Selection of balloon skin material 
<*-Shape and dimensions of balloon 

o Kinematics of the balloon system 
Altitude as a function or time 
Velocity as a function of time 
Acceleration as a function of time 
Balloon cycle design (Cooling instruments and 
transmission of data to eartn) 


o Dynamics of the balloon system 

Buoyancy force as a function of time and altitude 
Drag force as a function of time and altitude 
Volume as function of time and altitude 
Stability 
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o Valving 

(^Condensation rate of R30 in the secondary 
balloon J 

Gas valving rate 

o Heat exchange phenomenon of balloon 

K Computer simulation c>f flow conditions and 
heat transfer as a function of time using Fluent 
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a Convectional heat transfer 
b Conductive heat transfer 
c Radiative heat transfer 


Bulk temperature as a function of time and 
altitude 

^Overall Buoyancy as a function of bulk temperature 




Report Format 


This report includes these four main sections: 


f\. Technical Approach: 


o A complete discussion of the trajectory model, 
o A complete discussion of the heat exchange model. 

2. Difficulties involved and assumptions made: a 
discussion on how some hard-to-find data and equations were 
quantified and estimated. 
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3. Results: a presentation and discussion of the overall 
J performance and deficiencies of the balloon system. Errors 
mvolved are analyzed and discussed. 


,<f 


4. Conclusion and recommendation: 

~A summary of the possible effects of the outcome of the 
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rehrfement >n; 3 recommenc * atlon f° r further study and design 
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TECHNICAL APPROCH 


A Word About The Weather 

The Venus Atmosphere Model (0 to 70 km) 


The mostly carbon dioxide atmosphere (98%) acts as a huge 
thermal blanket, trapping most of the impinging thermal radiation 
from the sun Surface temperature is as high as 460 Celsius. Thick 
layers of sulphur dioxide cloud exist up to an altitude of about 60 km 
(Hunter 1983). Iimnediately above the cloud layer, the temperature 
drops to about -10 Celsius. Pressure and density also drops 
significantly (see Fig 8c\l - 8c.3 orythe page). 

Double Balloon System and Reversible Fluid 

In order to ensure a long duration operation, good altitude 
control design is essential . The life-span or the entire system 
depends on now well the onboard instruments are insulated and 
maintained. Temperature control of the interior of the craft is 
directly related to altitude control (buoyancy control). Our design 
uses a double balloon system. The primary balloon is filled witn , 

hy^Pgenwhile the secondary balloon uses a reversible refrigeranK_^A2f>' 
A< j- fig- 8b shows an estimation of the condensation rate of R Versus 
amtude. As reported by JPL^RJbegins to condense at an altitude of 
about 56 km above the surfacepf Venus. 


~ ,1 heijaUopn system continues to rise until thei^net buoyancy is 

less man its entire. weight. The liquid refrigerant is trapped in a 
small heat exchanger at the base or the secondary balloon. A closed 
valve prevents the R from evaporating back into the secondary 
balloon. The balloon system continues to descend below. tl^56krh — 
altitude until the valve is re-opened; the superheated ^evaporates 
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data to Earth is completed. 
Specifications 
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As suggested bv TPL, the balloon skin is to be manufactured out 
ot polyethyene. The balloon is a free pressure system. It is entirely 
h“ e e to expand and there are no stresses on its skin. J 

i-ui 00nS i are designed to haye a fixed diameter of four meters; 
only the lengths of the balloons change as the system ascends or 
descends. Its cigar shape is desienedfor drag reduction. See 
appendix M for detailed specification and dimensions of the 
balloon system. 
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CONDENSATION OF R30 IN 2ND BALLOON 

(ESTIMATION) 
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The Mathematical Model 


Earlier balloon users sought mathematically simple models for 
predicting balloon ascent rate. The earlier equations were purely 
dynamics equation. These equations did not include thermodynamic 
and heat transfer effects and were thus not very accurate. 

U 11 In Model, heat transfer effects within the balloon, on the 
balloon skin and outside the balloon,are carefully analyzed. Due to 
the size of the two balloons, the temperature of the entire system can 
neither be considered uniform nor equal to that of its immediate 
environment. A separate section is devoted entirely to the 
investigation of the transient heat transfer phenomenon of the 
balloon (see next section). The purpose of this investigation is to 
obtain a bulk interior temperature as a function of time. This bulk 
temperature is then mcottperated into the set of differential 
equations of motion to solve for its flight path. 

If we assume that the balloon trajectory problem is two- 
dimensional in nature, then the following differential equations 
define the model when the balloon is partially full; either when it is 
floating or when it is moving vertically, upward or downward. As 
mentioned earlier, the balloon is a free pressure system. Its bulk 
pressure is assumed to be that of its surroundings at all times. 

The differential equations: 

- (1) dv/dt + rhogVo11 

(2) V =dh/dt 

(3) Vol2 = nl'R'TbuUil/p (The gas law is assumed) 
where: 

V = velocity of balloon 
mtot = total mass of balloon system 
rhog = density of atmosphere at altitude h 
g = acceleration due to gravity 
Voll = volume of hydrogen balloon 
Vol2 = volume of K balloon 

SWK = bulk temperature of interior of hydrogen balloon 
Tbulk2 = bulk temperature of interior of K balloon 
I = external pressure at altitude h 
nhn2 = number of moles of gas 
Cd = coefficient of drag 
Ap = projected area ofballoon 
h = altitude of balloon 
R = gas constant 
k = Virtual mass coefficient 



The Drag Coefficient 


As advised by Jay Wu of JPL, the drag coefficient of the balloon was 
assumed to be O.b at all altitudes. _ 7 

^ . f t0 <> ) ■ 

The Added Mass (Virtual Mass) 


The conditions under which the added mass term applies are not too 
well defined. Generally, the added mass term becomes important 
when the density of the working fluid (inside the balloon) is lower 
than that of its environment’s (Poter 1993, P97). No data on the 
mass coefficient was available for balloon shape. We estimated the 
coefficient to tjaOOl^ ^ 

Solving the Differential Equation on Quattro Pro £ 

Equations (1), (2) and (3) were solved simultaneously using the 
Runge-Kutta Method of order 4. Equation (1) is a second order, non- 
linear, non-homogeneous differential equation. A step-wise#" 
iterative process was used to solve for the velocity ana displacement 
of the balloon system as a function of time. 

The time step was staggered from 0.05 sec to 4 sec in successive" 
blocks to ensure convergence. The initial conditions were: 

At time t = 0 ; 

displacement (t=0) = 40000m (balloon is released from 

entry vehicle) 

velocity (t = 0) = -11m/ s (as determined from entry 

* calculations) 

For a more detairaescription, see spreadsheet on appendix N - P. 
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HEAT TRANSFER 


Computer Simulation 


The commercial fluid dynamics software. Fluent was used to 
investigate the heat transfer that occurs between: 


a The balloon skin and the exterior of the balloon 

b) The balloon skin and the interior of the balloon 

c) The overall heat transfer; the bulk temperature 
^4, the balloon system as a function of time 


Because of the numerous variables involved and the 
complexity of the problem, a simplification was made. The balloon 
Was jreated as a 2-dim e nsional object. The problem's complexity 
arises from the tact 1) that the environmental conditions are J 

al^chan in a ^th^' an< ^ ^ *h e conditions within the balloons are 



— 0 _._ time. To further simplify the 

environmental temperature was held constantrTjnly u. 
conditions of the balloon system were permitted to vary 


i . A Grid of the balloon was initially created (see Fig 8fLThe 
inlet conditions were taken to be those found at an altitude ofS^km 
abpve the Venusian surface. This altitudeis crucial because this is the 
point where the R30 begins to condense.^The condensation rate will 
provide an msight mto now much and how soon the balloon loses 
buoyancy once it passes the {56 km point. 


/■m. T . Fluent uses the conventional fluid mechanics equations 

(J\lavier : Stokes) and heat transfer equations (Free and Force 
convection) to generate results. A radiation file was also created in 
Fluent to model the radiative effects of the balloon and its -e? 

surroundings. Heat transfer data foii^lJLlwasiSro\uaedT5vTavV^iT^ 
The overallTieat transfer coefficient of the balloon was found to be 
4.4 (For computation, see appendix Q ). Laminar flow was assumed 
because of the low velocity of the flow^/ 

Snapshots of the bulk balloon temperatures were taken 
every 15 minutes (see figures 8 g.l-.5). / * 

The whole program look about 8 hours to run and many more hours 
to create, ’ 

a i.- Fluent uses a finite difference method to compute results. 

A time step of 0.1 second was specified initially. 


s<rs - /2 3& 


the skin 


9H e a limitation of the resolution capabilities in Fluent, 
of the balloon cannot be made as thin as that specified in the 

6S12T1 snpri fi ra f i nn q Tn mmnonc^fn 
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actual design specifications. To compensate for this, the 
thermal conductivity was proportionally increased. / 
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Results 


The Table below shows the bulk temperature of the system as a 
function of time. J 


TIME 

(mins) 

HYDROGEN BALLOON 
(K) 

R30 1 

0 

315 

328 

15 

304 

314 

30 

300 

308 

60 

291 

301 

180 

291 

292 


(K) 


Table 4. The bulk temperature as a function of time. 

Fig 8h on the following page shows a plot of the above results. 
As expected, the R30 balloon took a longer period of time to 
equilibrate thermally (compared to the H2 balloon). 

Looking at the curves of the temperature vs. time plots, it can be 
mferred that the bulk temperature of the hydrogen balloon is about 
2Q3Lelsius higher than the surrounding during its asceh<3> and about 
20 Celsius lower during its descend ^ 

it % 

This is how it is approximated: 


At altitude = 56 km and at time t = 0 

Environment Temp. = 289 

H2 Balloon Temp, (initial guess) = 315 K 

Temperature difference =^5 degrees 

At time = 20 minutes later, balloon temperature has 

dropped to 303 Celsius. 

At altitude = 59 km and at time t = 20 minutes 
Based on an ascend rate of 3 m/s at that height) 
environment Temp. = 282 degrees 
H2 Balloon Temperature = 303 Celsius 
Temperature difference = 21 degrees 


At altitude = 64 km and at time t = 40 minutes 
(Based on an ascend rate of 4 m/s at that height) 
Environment Temp. = 276 degrees 
H2 Balloon Temperature = 296 degrees 
Temperature difference = 20 degrees 



IPERA' 









Hence it can be seen from the reasoning in the preceding pages 
that the difference between the interior bulk temperature of tne Pi2 

_ * • i 1 rs r\ 


balloon and the environment will converge to approximately 20 
degrees. Following the same line of reasoning, ft was found that the 
bulk temperature or the R30 balloon is about 22 degrees higher than 
the surrounding at all altitudes during the ascend cycle. 




} 






On the descend 
With this knowledge, 
spreadsheet (see app 
for the 'bulk H2 bailboj 
temperature'. The Rui 
durmg its iterative pr 
results into Runge-Kiitta process 



3 

cle, the temperature difference is reversed. rt 
ent back to the initial 
N) to create the 'columns' 
temperature' and 'bulk R30 balloon 
ge-Kutta algorit hpi a ddresses these cells 
cess. This is nov^wOnc^dperate the Fluent 

r 








Difficulties involved and assumptions made: 







The difficulties involved in modelling the balloon trajectory 
are numerous: J 


o Lack of data on condensation of R30 above 56 km 
o Lack of data on important heat transfer coefficients 
o Lack of precise formulas for calculating drag coefficients r 

o Lack of data on R30 (JPL was unable to forward any data ^ c 

on R30 to us; however, data on Freon R1 14 was forwarded ) c - u ' < ~ * 

instead.) v ‘ 

o Lack of references on space balloon trajectories 
o Lack of data on Venusian wind conditions 
o Scope of problem was too wide 
o Initial velocities had to estimated for the 
calculation of the Reynold's number 
o Lack of information on dimensions of system 
o Difficulties in solving a system of non-linear differential 
equations 

o Radiative properties of R30 were not available 
o No data on solar intensity below the 60 km was 
available 


Assumptions 


o 

o 

o 

o 

o 


o 


o 

o 

o 


Initial mass of various components were estimated 
3-dimensional problem was reduced to 2-dimensional 
for ease of treannent 
Wind forces were ignored 
Condensation rate of R30 was estimated 
Initial gas requirements (no. of moles required) 
were estimated 

Initial velocities were estimated 

Both internal and external flow were estimated 

to be laminar 

Balloon skin thickness was estimated 

Radiative properties of R30 were estimated 

Fluent has a great deficiency; it cannot 

model a balloon that changes in size constantly. ± -L ' 

As the balloon rises, its bulk temperature will lyA . ' 

decrease. If the volume of the balloon remains 

constant, the interior of the balloon will develop a negative 

pressure. In our model, this effect is 

ignored 
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RESULTS 


. The velocity of the balloon after the initial release stabilizes 
after 3 seconds (see Fie 8d.l). The balloon continues to rise up to a 
maximum height of 64 km. At this height it is far enough above the 
cloud layer to permit clear transmission (see Fig 8d.2) This occurs 
three hours after deployment. The maximum velocity attained bv 
the balloon was found to be 4 m/s Jsee Fig 8d.3). y 

., , When most of the gas in the secondary balloon has condensed, 

the downward forces predominate and the balloon system 

descends. The condensed fluid seeps into a heat exchanger at the 

base. As the balloon passes the 56 km altitude, th<£R^aporates / 
agam due to the increase in temperature. To prevent this, a valye^t 
the neck of the heat exchanger closes, trapping the condensect^I) 
ihe balloon continues to descend until the valve opens, releasing the 
8 e ^ e 3 ea t e d^< ^ or a more detailed trajectory description, see Figs 

V? 

The condensation rate of R30 had to be estimated. No data on 
this phenomenon was available. This had led to some discrepancies 
between the predicted result and the actual result generated bv the 
Runge-Kutta algorithm. Figure 8d.2 shows that the balloon 
remains afloatat an altitude of about 54 km after its initial 
deployment. The balloon is expected to descend as it loses buoyancv 

1 his discrepancy cannot be resolved until y ' 

more information on the condensation rate is made available to us. 



How long are ^^gom^tostay up there? Period of each 
cycle? At these point, these^nswerscannot be answered because JPL 
has not providecHUVith any specific requirements. The cycle 
period can be adjusted according to cooling needs of the 
instruments within the gondola and transmission requirements. 

At this point ( fCv?3 Ao not have any information on these requirements. 
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The Need for a Heat Exchanger 


u 




To accurately control the height of the balloon versus time a 
heat exchanger was added to the system. When the balloon reaches 
a height of 56 km, the gaseous methylene chloride, which was 
providing buoyancy, will condense into a liquid. This reduces the lift 
force and the balloon will fall. It is the job of the heat exchanger to 
capture the condensed liquid and provide enough heat where the 
liquid will turn back into a gas. A valve will control the flow of tha A 
vaporized or condensed liquid. yi&eJl / 

Vaporizing the Phase Changing Material &tj***~* y *'^ 

The first step is to calculate the energy required to vaporize the 
liquid. To find the heat of vaporization (AH V ) we can use the 2 point 
form Clausius - Clapeyron Equation. 

log (P 2 / Pi) = (AH V / 2.303 * R) * (1/Ti - I/T 2 ) 
where P is the pressure, 

T is the temperature, 

R is the gas constant. 

The 1 and 2 represent different states as found on the Van't Hoff 
plot. The heat of vaporization is found to be 31682.5 Joules / mole 

(Appendix R- U ). To vaporize all 200 mol«of methlyene chloride it 
would take 6336.5 KJ. 

Now(wB>vant to calculate the energy required to heat the fluid 
from 333K to 366 K. This energy is given by ^ ^ 

E = AT * heat capacity 10* ' 

The heat capacity for methylene chloride is 99997.6 J/degree^The 
total energy to heat the fluid is 3299.9KJ. Now^need to put in 
9636.4 KJ of energy it to a gas^ 

Size of Heat Exchanger and Stresses Involved 

One limiting factor in the size of the heat exchanger is the 
space that it will require inside of the entry vehicle. For this reason it 
must be smaller than 13.8 inches long . By knowing the density of the 
liquid and the mass of theTlquid, the volume can be found. In this 
case it was 0.013 m 3 . N/xt an optimizing sequence was done to find 
the best radius, length, And thickness that would hold the liquid and 
resist the pressure involved from the phase change (Appendix V-X ). 

" 0 
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In this analysis of the stresses we will treat it as a thin walled 
cylindrical pressure vessel. The stress that we are concerned with is 
the hoop stress given by 

< n = F * r/t sfr2 ^ 

Where P is the gage pressure, Z 

r is the radius, j 

t is the thickness. 

The gage pressure is the difference in pressure between the inner 
pressure and the atmospheric pressure. For this model we will be 
using a n aluminum shel l with a yield stress of 416 M Pa. The final 
dimensions of the heat exchanger is a 0.1143 m inner radius cylinder 
0.3117 m in length with a thickness of 0.0138 m. The final mass of this 
unit is 8.5 kg. 

t/ 

Concept of a Heat Exchanger 

At low altitudes we have a hot fluid, the carbon dioxide 
atmosphere flowing across a vertical tube which contains the phase 
changing material, methylene chloride. The temperature difference 
between the atmosphere and the inner fluid will cause energy to 
transfer and heat up the methylene chloride changing it from a liquid 
to a vapor. The three types of heat transfer present are conduction, 
radiation, and convection. For all three types of transfer we 
assumed a one dimensional , steady state , vertical cylinder without 
fins. Other assumptions made are noted in the sections. 

Conduction 

Conduction will tell how much heat is lost through the walls of 
the heat exchanger. 

This is given by 

q = -k A (dT / dx) 
where q is the heat rate in Watts, 
k is the thermal conductivity of the material, 

A is the exposed area, 

dT/dx is the temperature gradient (Incropera). 

The negative sign shows that heat is transferred to the cooler side of 
the material. For a thin member we can assume this to be negligible 
but in our case this was calculated. In reality this is a two 
dimensional transient problem but due to the slow velocity, the 





material will have time to come to equilibrium with the 
surroundings. The calculations to^rovelhis assumption can be 
found in Appendix Y. ^ 

JessriFf' sw'T- y*&<~f***) 


Radiation 
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Radiation will tell how much heat is gained and reflected to the 
atmosphere. The general equation is given by 
q = A e a (T s 4 - T sur 4) 
where q is the heat rate in Watts, 

A is the exposed area, 
e is the emissivity of the material, 
a is the Stefan-Boltzmann constant, 

T s and T sur are the temperatures of the 
atmosphere and surface of the material 
respectively (Incropera). 

This equation will give us the heat rate from the atmosphere to the 
balloon, but it can also give us the heat rate from the balloon to the 
methylene chloride gas by changing the subscripts on the 

temperature to surface and gas respectively. We also need to ^ , 0 / 

consider the amount that is absorbed from solar radiation. The solar 
flux Venus receives is 2600 W/m 2 of which 132 W/m 2 is absorbed. It 
is the latter number that we use to calculate how much solar flux is 
absorbed by the heat exchanger. The calculations can be found in 
Appendix Z. 


Forced Convection 


This is the greatest source of heat transfer to the cylinder. The 
heat rate for convection is given as ^ 

q = h A (T sur - T s ) 

where q is the heat rate in Watts, t/^ r ^ 


/>?( 


A is the exposed area, _ 
h is the convection coefficient, 

T s and T sur are the temperatures of the 


>tf£> 




atmosphere and surface of the material 
respectively (Incropera). 





In order to find the convection coefficient we need to calculate the 
Reynolds (Re) and Nusselt (N u ) numbers. This will be treated as 
external flow over a cylinder where the following equations apply 
Re = V D / v 
N u = C Re m Pr(l/3) 
h = N u k/D 

where V is the velocity of the wind, 

D is the diameter of the tube, 
v is the viscosity of the atmosphere, 

Pr is the Prandtl number, 
k is the thermal conductivity. 

C and m are constants dependent on the Reynolds number. 
Numerical calculations can be found in Appendix AA - AB. 




Free Convection 

This is a source of heat transfer/to the inside gas. The same 
general equation applies to free convection. The difference is in the 
method to find the convection coefficient, h . This is modeled as 
laminar flow on a vertical isothermal plate where the following 


r 


equations apply 

Gr = g P (T sur - Tg as ) / v2 
Nu = (4g(Pr)/3)(Gr/4)(V4) 
f h = Nu k/D 

where Gr is the Grashof number. 
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is the volume coefficient of expansion (Incropera). mu* 

fnr tVlic ni 1 anf iH r ^1. ~ £ ' 
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Calculationsjo r this quantity were not d one due to the lack of ' s 0 ”* 

i nformation of the properties of methyle n e chloride . Inside the heat 




exchanger the liquid is boiling and thus the above equations are no 
longer valid. The equation for pool boiling is 

q" = mi hfg [ g (pi - pv) / o ] 1/3 [Cpl ATe/C S/ f (hf g ) P r n] 3 J 
where q" is the heat flux, 
ml is the mass of the liquid, 
hfg is the latent heat of vaporization, 
o is the surface tension, 

L p is the density. 

The letters v and 1 represent the vapor and liquid states, C S/ f and n 
are constants (Incropera). (See Appendix AC) 
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Problems With Heat Transfer 


By doing all of the above calculations, the balloon will receive 
enough energy to heat and vaporize the liquid that has condensed 
after two hours. All the calculations done have been based on a 
steady state model floating between the altitude of 48 to 60 km. This 
time would actually be longer due to the fact that this is a transient 
problem and the heat rate is proportional the temperature ^ 
difference. Another problem is calculating the horizontal winds. The 
graph in Appendix shows a broad spectrum for the velocity of the 
wind. The stress calculations show that the heat exchanger will be 
safe for cycling about the 56 km mark but might not withstand the 
pressure near the surface. Future work should be done to £nHy~ 
calculate the stresses depending on what range the balloon system 
will see.' / The material for the heat exchanger must be lightweight, 
have a high thermal conductivity and also be able to withstand the 
pressure near the surface. „ / \ 



CONCLUSION 


(fWihelieve thatjomjnvestigation and design as it is defined 
will provide an added source of literature to the Venus balloon 
concept. Our preliminary study will lay the ground work to further 
refinement of the trajectory model. W e Tecommend that 
actual heat transfer experiment be conducted and results compared 
against those generated by computer simulation. 
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Department 

Professor G. Meyers, at UW-Madison, Mechanical Engineering 
Department ( Thermo ) 

Professor Mitchell, at UW-Madison, Mechanical Engineering 
Department ( Heat Transfer) 
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Department ( Fluid ) 
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Astronautics Department 
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The Venus Balloon used the following computer software packag es 
for modeling . 

Excel 

Cricket Graph 
Quattro Pro 
Fluent 

Mathematica 
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Appendix A: Definition of Symbols 


m 


A 


Bo 

Cd 

Cf 

Dh 

Ein 


Etot 

e 

e 

g 

Ho 

hvap 

Kc 

Kd 

K g 

kchar 

ksheU 

kvirgin 

L 

M 

m 

Nu 

oyield 

P 

Pr 

V 

pchar 

^virgin 

Poe 


Exposed area (length^), assumed constant 
Aerodynamic deceleration (length /time^) 

Energy parameter (mass*length2/time2) 

Drag coefficient (unitless), assumed constant 
Skin Friction coefficient (unitless) 

Characteristic length (length) 

Convective energy seen by the craft 
(mass*length2 / time2) 

Total kinetic energy lost by the craft 
(mass* length2 / time2) 

Euler constant, 0.5772156 (unitless) 

Thermal strain (unitless) 

Gravitational acceleration (length/time^) 

Scale height (length) „ — ^ ' 

Heat of shield vaporization (length^/time^) 

Ratio of specific heats (unitless), assumed constant 
Drag parameter (length 3 /mass) 

Universal gas constant divided by molecular weight 
(Iength2*temperature/time2) 

Thermal conductivity of the char layer (length) 

Thermal conductivity of the aeroshell (length) 

Thermal conductivity of the shield (length) 

Vehicle length (length) 

Mach number (unitless) 

Spacecraft mass (mass), assumed constant 
Nusselt number (unitless) 

Yield stress (mass / length-time^) 

Dynamic pressure (mass / length-time^) 

Prandtl number (unitless) 

Atmospheric density (mass /length 3 ) 

Char density (mass /length 3 ) 

Shield density (mass /length 3 ) 

Atmospheric density at stagnation point 
(mass /length 3 ) 

Convective heat flux (mass /time 3 ) 





R 

Universal gas constant (length 2 /time 2 -temp) 

Rn 

Nose radius of curvature (length) 

Re 

Reynolds number (unitless) 

r 

Base radius (length) 

T 

Atmospheric temperature (temperature) 

t 

Aeroshell thickness (length) 

tchar 

Thickness of the shield char layer (length) 

tshell 

Thickness of the aeroshell (length) 

tvirgin 

Thickness of the untouched shield (length) 


Viscosity (mass /length- time) 

V 

Current velocity (length /time) 

Vchar 

Volume of the char layer (length^) 

Vshell 

Volume of the aeroshell (length^) 

V virgin 

Volume of the shield (length^) 

Vo 

Entry velocity, 11540 meters per second 



m 


b v 



Appendix B: Specific Model of a Venusian Atmosphere 

Most models of planetary atmospheres use an isothermal 
constant composition algorithm. This facilitates the creation of a 
fictitious reference height that makes the calculations of 
acceleration, velocity, and position very easy. A summary of the 
relevant equations for this case is shown below (Weisel). An 
improvement on the model is to remove the isothermal constraint, 
adding a lapse rate term to the equations of motion. The lapse rate is 
the change in temperature with respect to the change in height. To 
keep calculations manageable, the atmosphere is usually divided into 
sections that have constant lapse rates (Regan). The designer has not 
settled for either of these approximations. The model upon which the 
design is based uses given values of pressure, temperature, and 
density at each altitude, plus correlations for gravity, viscosity, and 
chemical composition. A section of the spreadsheet model is shown in 
Figure B. From this data, a more realistic profile for deceleration can 
be created. Figures Ba-Bd show comparisons between typical plots 
from the spreadsheet and general plots of expected entry vehicle 
results. Accuracy could be further improved by allowing angle, mass, 
and drag coefficient to vary during entry instead of remaining 
constant. This model can provide accurate profiles for any ballistic 
entry to a reasonable degree. It is interesting to wonder how this 
model would fare against trajectory simulation software. 

SUMMARY OF DEVELOPED EQUATIONS 
Cd*A*Ho 

Kd = H 0 = — 

m*sin (angle) g 

B 0 = l/2*m*V 0 2 - m*g*H 0 *[e + ln(Kd*p)] 

l/2*m*v2 = [B 0 + m*g*H 0 *(e + ln(Kd */?))]* exp[-Kd*p] 

-Kd* sin (angle)*B 0 

a = * p * exp[ - Kd*p] 


m*H 0 


6^1 



o 


Height (km) 

Height (m) 

Temperature 

Bluntcricket ^ 

7 fa 

Gravity i Density 

j s+C 

L Angle 

/ 37 

-1.274 

1 180.000 

1 80000.000 

125.000 

J ! ■ 

8.365 

0.000 

2 176.000 

1 76000.000 

125.000 

8.376 

0.000 

-1.274 

3 172.000 

1 72000.000 

125.000 

8.387 

0.000 

-1.274 

4 168.000 

168000.000 

126.000 

8.397 

0.000 

-1.274 

5 164.000 

164000.000 

126.000 

8.408 

0.000 

-1.274 

6 160.000 

1 60000.000 

127.000 

8.419 

0.000 

-1.274 

7 156.000 

1 56000.000 

1 29.000 

8.430 

0.000 

-1.274 

8 152.000 

1 52000.000 

131.000 

8.441 

0.000 

-1.274 

9 148.000 

148000.000 

134.000 

8.452 

0.000 

-1.274 

1 0 144.000 

144000.000 

140.000 

8.462 

0.000 

-1.274 

1 1 140.000 

140000.000 

147.000 

8.473 

0.000 

-1.274 

1 2 136.000 

136000.000 

146.000 

8.484 

0.000 

-1.274 

13 132.000 

1 32000.000 

142.000 

8.495 

0.000 

-1.274 

1 4 1 28.000 

1 28000.000 

135.000 

8.506 

0.000 

-1.274 

1 5 124.000 

124000.000 

125.000 

8.517 

0.000 

-1.274 

16 120.000 

120000.000 

1 17.000 

8.528 

0.000 

-1.274 

1 7 1 16.000 

1 16000.000 

1 90.200 

8.540 

0.000 

-1.274 

1 8 114.000 

1 14000.000 

1 90.200 

8.545 

0.000 

-1.274 

1 9 1 12.000 

1 12000.000 

190.200 

8.551 

0.000 

-1.274 

20 1 10.000 

1 10000.000 

194.400 

8.556 

0.000 

-1 .274 

21 108.000 

1 08000.000 

196.300 

8.562 

0.000 

-1.274 

22 106.000 

1 06000.000 

190.100 

8.567 

0.000 

-1.274 

23 1 04.000 

1 04000.000 

1 83.000 

8.573 

0.000 

-1.274 

24 1 02.000 

102000.000 

1 80.700 

8.578 

0.000 

-1.274 

25 100.000 

100000.000 

174.100 

8.584 

0.000 

-1.274 

26 98.000 

98000.000 

173.300 

8.590 

0.000 

-1.274 

27 96.000 

96000.000 

1 68.900 

8.595 

0.000 

-1.274 

28 94.000 

94000.000 

168.900 

8.601 

0.000 

-1.274 

29 92.000 

92000.000 

167.100 

8.606 

0.001 

-1.274 

30 90.000 

90000.000 

165.400 

8.612 

0.001 

-1.274 

31 88.000 

88000.000 

165.200 

8.618 

0.002 

-1.274 

32 86.000 

86000.000 

163.800 

8.623 

0.003 

-1.274 

33 84.000 

84000.000 

164.700 

8.629 

0.005 

-1.274 

34 82.000 

82000.000 

1 81 .1 00 

8.634 

0.008 

-1.274 

35 80.000 

80000.000 

195.800 

8.640 

0.012 

-1 .274 

36 79.000 

79000.000 

204.100 

8.643 

0.014 

-1.274 

37 78.000 

78000.000 

212.400 

8.646 

0.017 

-1.274 

38 77.000 

77000.000 

215.100 

8.649 

0.021 

-1.274 

39 76.000 

76000.000 

216.900 

8.651 

0.025 

-1.274 

40 75.000 

75000.000 

221 .100 

8.654 

0.031 

-1 .274 

41 74.000 

74000.000 

227.300 

8.657 

0.036 

-1.274 

42 73.000 

73000.000 

227.800 

8.660 

0.044 

-1 .274 

43 72.000 

72000.000 

233.500 

8.663 

0.053 

-1.274 

44 71.000 

71000.000 

234.900 

8.665 

0.063 

-1.274 

45 70.000 

70000.000 

243.500 

8.668 

0.074 

-1.274 

46 69.000 

69000.000 

230.500 

8.671 

0.094 

-1.274 

47 64.830 

64830.000 

246.300 

8.683 

0.1 96 

-1.274 

48 64.000 

64000.000 

247.000 

8.685 

0.228 

-1.274 

49 63.000 

63000.000 

249.300 

8.688 

0.272 

-1.274 

50 62.000 

62000.000 

254.800 

8.691 

0.318 

-1.274 

51 61.000 

61 000.000 

259.700 

8.694 

0.373 

-1.274 

52 60.000 

60000.000 

268.100 

8.697 

0.429 

-1.274 

53 59.000 

59000.000 

272.000 

8.700 

0.500 

-1.274 

54 58.000 

58000.000 

276.000 

8.702 

0.583 

-1.274 

55 57.000 

57000.000 

281 .300 

8.705 

0.673 

-1.274 

56 56.000 

56000.000 

288.900 

8.708 

0.769 

-1.274 


Sun, May 8, 1994 2:10 PM 
(rads) Mass 


1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 

150.000 
1 50.000 

150.000 
1 50.000 
1 50.000 
1 50.000 

150.000 
1 50. 

150. 

1 50/000 
1 5CX000 
1 5 
1 50 
1 50. 

1 50. od 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 

150.000 
1 50.000 

150.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 
1 50.000 


(kg) 

*7 


& 


n l" 





9 ** 






,lw * '\ Bluntcricket 

Sun , 

May 8 , 1994 2:10 PM 

Area 

( m A 2 ) Cd 

k * T/g 

Ho 

Kd 

£olumn 1_3 

•) Deceleration 

1 1.214 

1 .062 

277.143 

4141 .407 

37.220 

1 1540.000 

0.000 

2 1.214 

1.062 

272.600 

4068.287 

36.563 

1 1540.000 

0.000 

3 1.214 

1 .062 

268.203 

3997.530 

35.927 

1 1540.000 

0.000 

4 1 .214 

1 .062 

259.822 

3898.572 

35.037 

1 1540.000 

0.000 

5 1.214 

1 .062 

251.948 

3775.573 

33.932 

1 1540.000 

0.000 

6 1.214 

1.062 

244.538 

3688.858 

33.153 

1 1540.000 

0.000 

7 1.214 

1 .062 

237.551 

3635.208 

32.671 

1 1540.000 

0.000 

8 1.214 

1 .062 

230.953 

3584.401 

32.214 

1 1540.000 

0.000 

9 1.214 

1.062 

224.71 1 

3562.794 

32.020 

1 1540.000 

0.000 

1 0 1 .214 

1.062 

218.797 

361 9.691 

32.531 

1 1540.000 

0.000 

11 1 .214 

1 .062 

213.187 

3698.442 

33.239 

1 1540.000 

0.000 

12 1 .214 

1 .062 

207.857 

3576.825 

32.146 

1 1540.000 

0.001 

13 1.214 

1 .062 

202.788 

3389.594 

30.463 

1 1540.000 

0.002 

14 1.214 

1 .062 

199.384 

3164.309 

28.438 

1 1540.000 

0.007 

15 1 .214 

1 .062 

1 97.021 

2891 .455 

25.986 

11540.000 

0.029 

1 6 1 .214 

1.062 

195.631 

2683.818 

24.120 

1 1540.000 

0.130 

17 1.214 

1 .062 

1 95.631 

4357.271 

39.160 

1 1540.000 

0.81 1 

18 1 .214 

1 .062 

194.259 

4323.924 

38.860 

1 1540.000 

1 .288 

1 9 1 .214 

1 .062 

194.259 

4321 .1 19 

38.835 

1 1540.000 

2.058 

20 1 .214 

1.062 

1 93.356 

4393.144 

39.482 

1 1540.000 

3.172 

21 1.214 

1 .062 

1 93.356 

4433.202 

39.842 

1 1540.000 

4.958 

22 1 .214 

1 .062 

193.356 

4290.395 

38.559 

1 1540.000 

8.096 

23 1 .214 

1 .062 

193.356 

4127.472 

37.095 

1 1540.000 

13.540 

24 1 .214 

1 .062 

193.356 

4072.948 

36.605 

1 1540.000 

22.414 

25 1 .214 

1 .062 

192.461 

3903.480 

35.082 

1 1540.000 

38.472 

26 1.214 

1 .062 

192.461 

3883.017 

34.898 

1 1540.000 

64.573 

27 1 .214 

1 .062 

192.461 

3781 .968 

33.989 

1 1540.000 

1 11.662 

28 1 .214 

1 .062 

192.461 

3779.508 

33.967 

1 1540.000 

1 88.783 

29 1.214 

1.062 

192.461 

3736.796 

33.583 

1 1540.000 

324.392 

30 1 .214 

1 .062 

191.574 

3679.338 

33.067 

1 1540.000 

554.407 

31 1.214 

1.062 

191.574 

3672.496 

33.006 

1 1540.000 

935.084 

32 1.214 

1 .062 

191.574 

3639.001 

32.705 

1 1540.000 

1 556.255 

33 1.214 

1 .062 

191.574 

3656.61 1 

32.863 

1 1540.000 

2503.387 

34 1 .214 

1 .062 

191.574 

401 8.098 

36.1 12 

1 1540.000 

3415.717 

35 1.214 

1.062 

191.574 

4341.417 

39.017 

1 1540.000 

4277.213 

36 1.214 

1 .062 

191.397 

451 9.809 

40.621 

1 1 540.000 

4575.1 62 

37 1 .214 

1.062 

191.397 

4702.079 

42.259 

1 1540.000 

4745.451 

38 1.214 

1.062 

191.397 

4760.297 

42.782 

1 1540.000 

4886.656 

39 1 .214 

1 .062 

191.397 

4798.566 

43.126 

1 1 540.000 

4858.835 

40 1 .214 

1.062 

191.397 

4889.888 

43.947 

1 1 540.000 

4558.937 

#* 41 1.214 

1 .062 

191.397 

5025.368 

45.164 

1 1540.000 

4016.905 

42 1 .214 

1 .062 

191.397 

5034.778 

45.249 

1 1540.000 

3408.794 

43 1 .214 

1.062 

1 91 .397 

5159.073 

46.366 

1 1540.000 

2627.004 

44 1 .214 

1 .062 

191.397 

51 88.310 

46.629 

1 1540.000 

1 885.626 

45 1 .214 

1 .062 

191.397 

5376.505 

48.320 

1 1540.000 

1191.921 

46 1 .214 

1 .062 

191.397 

5087.801 

45.725 

1 1540.000 

719.051 

47 1.214 

1 .062 

191.397 

5429.148 

48.793 

1 1540.000 

7.728 

48 1 .214 

1.062 

191.397 

5443.100 

48.918 

1 1540.000 

1 .849 

49 1 .214 

1.062 

191.397 

5491.988 

49.358 

1 1540.000 

0.234 

50 1 .214 

1.062 

191.397 

561 1.316 

50.430 

1 1540.000 

0.019 

51 1.214 

1 .062 

191.397 

5717.355 

51 .383 

1 1540.000 

0.001 

52 1 .214 

1 .062 

191.397 

5900.352 

53.028 

1 1540.000 

0.000 

53 1 .214 

1 .062 

191.397 

5984.225 

53.782 

1 1 540.000 

0.000 

54 1.214 

1 .062 

191.397 

6070.241 

54.555 

1 1540.000 

0.000 

55 1 .214 

1 .062 

191.397 

61 84.782 

55.584 

1 1540.000 

0.000 

5 6 1.214 

1 .062 

191.397 

6349.800 

3 - ( 

57.067 

1 1540.000 

0.000 



Bluntcricket 


Sun, May 8 , 1994 2:10 PM 


Deceleration 

1 0.000 
2 0.000 

3 0.000 

4 0.000 

5 0.000 

6 0.000 

7 0.000 

8 0.000 
9 0.000 

1 0 0.000 
1 1 0.000 
1 2 0.000 
1 3 0.000 
1 4 0.001 
1 5 0.003 
1 6 0.013 
1 7 0.083 
1 8 0.131 
1 9 0.210 

20 0.324 

21 0.506 

22 0.826 

23 1.382 

24 2.287 

25 3.926 

26 6.589 

27 1 1.394 

28 1 9.264 

29 33.101 

30 56.572 

31 95.417 

32 158.802 

33 255.448 

34 348.543 

35 436.450 

36 466.853 

37 484.230 

38 498.638 

39 495.799 

40 465.198 

41 409.888 

42 347.836 

43 268.062 

44 192.411 

45 121.625 

46 73.373 

47 0.789 

48 0.189 

49 0.024 

50 0.002 

51 0.000 

52 0.000 

53 0.000 

54 0.000 

55 0.000 

56 0.000 


7 ACWOU 

Bo j Column 17 J V A 2 V 

1 01 13000000 . 99878700 et) 0 133171600.00 1 1540.000 
10108000000 . 9987870000.0 133171600.00 1 1540.000 
10104000000 . 9987870000.0 133171600 . 00 1 1540.000 
1 0098000000 . 9987870000.0 1 331 71 600.00 1 1 540.000 
10093000000 . 9987870000.0 133171600.00 1 1540.000 
10088000000 . 9987870000.0 133171600.00 1 1540.000 
10084000000 . 9987870000.0 133171600.00 11540.000 
10080000000 . 9987870000.0 133171600.00 11540.000 
10076000000 . 9987870000.0 133171600.00 1 1540.000 

10074000000 . 9987870000.0 133171599.00 11540.000 
10071000000 . 9987870000.0 133171599.00 11540.000 

10064000000 . 9987870000.0 133171596.00 11540.000 

1 0055000000 . 9987870000.0 133171586. 00 11 539.999 


10046000000. 9987870000 

1 0036000000. 9987870000 

10027000000. 9987870000 

10039000000 . 9987870000 . 
10036000000 . 9987870000 . 

10034000000 . 9987870000 . 
10032000000 . 9987870000 . 
10030000000 . 9987870000 . 
10026000000 . 9987870000 . 
10022000000 . 9987870000 . 
10019000000 . 9987870000 . 
10015000000 . 9987870000 . 
10013000000 . 9987870000 . 

10009000000 . 9987870000 . 
10007000000 . 9987870000 . 
10004000000 . 9987870000 . 

10001000000 . 9987870000 . 

9998726352.0 9987870000 

9996085598.0 9987870000 

9993514464.0 9987870000 

9991360010.0 9987870000 

9988948470.0 9987870000 

9987678998.0 9987870000 

9986352625.0 9987870000 

9985034942.0 9987870000 

9983708216.0 9987870000 

998231 7404.0 9987870000 

9980846690.0 9987870000 

9979539623.0 9987870000 

9978029122.0 9987870000 . 

9976668536.0 9987870000 . 

9974945436.0 9987870000 . 

9974369662.0 9987870000 . 

9967807649.0 9987870000 . 

9966668249.0 9987870000 . 

9965160683.0 9987870000 . 

9963341776.0 9987870000 . 
99615521 12.0 9987870000 . 

9959375786.0 9987870000 . 

9957650679.0 9987870000 . 

9955889390.0 9987870000 . 

9953955595.0 9987870000 . 

9951713072.0 9987870000 . 


0 1 331 71 554.00 1 1 539.998 
0 133171427.00 11539.993 
0 133170874.00 11539.968 
0 133164247.00 1 1539.681 
0 133160008.00 1 1539.498 
0 133153087.00 11539.198 
0 1 33142579.00 1 1 538.743 
0 133125818.00 1 1538.016 
0 13309921 7.00 1 1536.863 
0 133055068.00 11534.950 
0 132981 136.00 1 1531 .745 
0 1 32858021 .00 1 1 526.405 
0 132647481.00 11517.269 
0 1 32287368.00 11501 .625 
0 131673766.00 1 1474.919 
0 130615885.00 1 1428.731 
0 128840012.00 1 1350.771 
.0 1 25788863.00 1 1215.563 
.0 120737994.00 10988.084 
.0 1 1231 1633.00 10597.718 
.0 99934499.000 9996.724 
.0 8377311 1 .400 91 52.765 
.0 74557896.200 8634.691 
.0 6481 5821 .600 8050.827 
.0 54765326.700 7400.360 
.0 44533465.400 6673.340 
.0 34643635.100 5885.884 
.0 25648876.900 5064.472 
.0 1 7900971.700 4230.954 
.0 1 1642350.200 3412.089 
0 6929935.480 2632.477 
03762719.200 1939.773 
0 1774182.890 1331.985 
0 9174.546 95.784 

0 1889.284 43.466 

0200.649 14.165 

0 14.228 3.772 

0 0.639 0.799 

0 0.017 0.132 

0 0.000 0.017 

0 0.000 0.001 

0 0.000 0.000 

0 0.000 0.000 


Time elapsed Total time 


0.347 

0.347 

0.347 

0.693 

0.347 

1 .040 

0.347 

1 .386 

0.347 

1 .733 

0.347 

2.080 

0.347 

2.426 

0.347 

2.773 

0.347 

3.120 

0.347 

3.466 

0.347 

3.813 

0.347 

4.159 

0.347 

4.506 

0.347 

4.853 

0.347 

5.1 99 

0.347 

5.546 

0.173 

5.719 

0.173 

5.893 

0.173 

6.066 

0.173 

6.239 

0.173 

6.413 

0.173 

6.586 

0.1 73 

6.759 

0.173 

6.933 

0.174 

7.106 

0.174 

7.280 

0.1 74 

7.454 

0.175 

7.629 

0.176 

7.805 

0.177 

7.982 

0.180 

8.1 62 

0.185 

8.347 

0.194 

8.541 

0.209 

8.750 

0.1 12 

8.863 

0.120 

8.983 

0.129 

9.112 

0.142 

9.254 

0.1 59 

9.413 

0.183 

9.596 

0.215 

9.81 1 

0.262 

10.073 

0.331 

10.404 

0.437 

10.841 

0.61 1 

11.453 

5.818 

17.271 

8.489 

25.760 

1 1.087 

36.847 

1 2.049 

48.896 

13.040 

61 .936 

14.051 

75.986 

15.119 

91.106 

16.318 

1 07.423 

17.572 

1 24.995 

18.837 

143.832 

20.066 

163.898 


^ X 



For each height listed in the spreadsheet, the values of Cd, A, 
m, angle, k, T, g, p, V 0 , and e are given. There are fifteen 
parameters, 10 knowns, 5 unknowns, and five equations. This set of 
equations has a closed form solution. The simple textbook model 
substitutes in for p: . : ’ ' 

p = p 0 * exp(-H/H 0 ) 

This is the exponential atmosphere equation and is only valid for 
isothermal constant composition atmospheres. 




Figure fta: Typical entry vehicle plots 



Velocity profile with re&pect to height for 
ballistic entry, t 4 degrees from horizontal 



entry 74 degrees from horizontal 



FiguresBb and 6c: Plots of designed vehicle parameters 
(Compare to a) 



convective flux profile for ballistic model 


with entry at 74 degr 


200 


from horizontal 
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Convective flux ( W/m^ 2 ) 


Figure 6d: Another plot to compare to(?a; FigsBb-W 
are similar to those shown in Figure Sa 


Appendix C: Numerical Justification of Shape and Material 
Selection 

i Pteflsc peruse t h e tables below for the four different shaped 
evaluated in thiG docign. Each table is based on calculations for a 250 
kg vehicle with the areas and drag coefficients for each shape lis ted 
It was desired to maximize the angle for which a 500g entry would 
occur, while simultaneously minimizing the peak heating rate, the 
total energy transfer, and the total time of entry. From these and 
other tables, the shape selection was narrowed down to the two 
highlighted choices. The worst heat shield material candidate 
properties were plugged into the conductive heat transfer equations 
to fmd which shape required less heat shielding. Some sample 
calculations follow, from which the reader should convince 
themselves that the best shape is the blunted cone. The second page 
of calculations documents the selection process of the shield and shell 
materials. Although there were only a small scope of materials 
considered, each one of these materials have been candidates for 
Venus missions before. 
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Appendix D: The Heating Problem 
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The design of the heat shield and aeroshell was influenced by / 
the thermal balance. Figurqy^hows the simplified version of the / 
problem. A model of one dimensional steady-state heat transfer W as 
used, incorporating the three knowns of desired interior 
temperature, maximum heating rate, and total convective energy ul ^ 

transferred. The surface temperature of the heat shield was yur* 

estimated based on the melting point of the material. The ratio of' / "~ i * * 
charred material to virgin material at any time was Jiaggledgyer, 
and node temperatures were set based on the established conditions. 

Pluggin g in all the above values gave the needed thickness of the 

at shield for thermal loading. This value was plugged in to volume / 

calculations to find the mass required. Finally, the calculated mass 
was checked against the total energy transfer term to see if all the 
material would be vaporized. A sample calculation followsy^ 
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Appendix E: The Loading Problem 


The design of the aeroshell was affected by the high loadings 
during descent. A modified approach based on the pressure vessel 
analogy was used to estimate a suitable thickness. The dynamic 
pressure was used as the value of the gage pressure in the 
calculations. A worst case scenario was developed and the thickness 
solved for. This thickness was compared to the thickness calculated 
for thermal loading, and the greater of the two values became the 
final result. A sample calculation follows: 
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Appendix ^ 


GxoftduoQ-c*. 



Equipment 


charge regulator 


power distributor 


transmitter 


command & data storag 


sensor * 


battery * 


camera * 


phase change material 


cable 1 


1 0 cable 2 


1 1 cable 3 


Dimensions L*W*H (in. A 3) 


5x4x3 


4 x 1.8x6 


4.7 x 7.1 x 8.72 


, 8 x 4.3 x 7 


2 


3.2 


2.5 


1.5 x 1.5 x 1.5 


5 m Length, 0. 01 5m radius 


5 m Length, 0. 01 5m radius 


5 m Length, 0. 01 5m radius 


^ 750 (J AS A? 3 




Mass (Kg) 


Power (W) 


8.0 




3 

1.5 

0 



Assupmtion: no information , The numbers were determined from the Viking lander 
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Power Generator ( Heat pipe) 


Foam 


Insulation 

section 


charge I 

regulator | 


ransmltter 



command 
& data 
storage 


material 


sensor battery camera 




Outer shell 


Fiber optic 


Front View 


Ayis~u57^ 



3 dimension 


2 dimension 
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MODEL OF THERMOCOUPLE DEVICE 





<{ V 







Appendix l— 


Hecct Prpe 





Container section 

Insulation section 
Pipe section 

Main container 
section 


L xW xH = 0.1 x 0.1 xO.1 = 1x10*3 mA3 

LxWxH = 0.00085 m** 

0.03 m -diameter , 

LxWxH = 0.1 5x 0.15 x 0. 32 = 7.2x1 OA3mA3 


Valve 



Thermo Couple 
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A PP^AJDIV X 


PHASE TRANSITIONS 


Clausius - Clapeyron Equation (2 point form) 

log (P 2 / Pi) = (AH vap / 2.303 R) [(1 /Ti)-(1/T 2 )] 

This equation relates temperature and pressure with the heat of 
vaporization. We can use values of pressure and temperature off the 
Van’t Hoff plot for a specific chemical. Units of AH va p — > kj / mol 

The heat needed for the vaporization of a liquid is called the heat of 
vaporization (or enthalpy of vaporization) 

HEAT CAPACITY 

Heat Capacity, C, is the quantity of heat needed to raise the 
temperature of the sample of the substance one degree C. 

q = C * AT 

where T = Tf - Ti and the units of C are Joule/ °C 


SPECIFIC HEAT 

The specific heat of a substance is the quantity of heat required 
to raise the temperature of one gram of that substance by one degree 
C. 

q = sp.ht. * mass * AT 

where T = Tf - T{ and the units of sp.ht. are Joule/ gram °C 


CRITICAL TEMPERATURE AND PRESSURE 

The temperature above which the liquid state of a substance no 
longer exists is called the critical temperature. The vapor pressure at 
the critical temperature is called the critical pressure. It is the 
minimum pressure that must be applied to a gas at the critical 
temperature to liquefy it 
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APPENDIX S 

Properties of Methylene Chloride (CH 2 CI 2 ) 

Taken from: CRC Handbook of Chemistry and Physics 67th Edition 


Pg C349 No. 9060 

Molecular weight 84.93 
Boiling Point 40 °C 

Melting Point -95.1 °C 
Density 1.3266 

Pg C671 Heat of Vaporization AHv = 7572.3 gram calories/gram 
mole (To convert to Joule/ gram mole, multiply listed by 
4.184) 

Pg C716 Limits of Superheat of pure liquids 
P = 0.101 MPa 
T = 394.8 K 

Pg E4 Thermal Conductivity k = 0.0002908 
[t(°C) = 0, t(°F) = 32] 

Pg F63 Critical Temperatures and Pressures 

Critical Temperature Tc (°C) 273 
Critical Pressure Pc (atm) 60 

Pg D59 Selected Values of Chemical Thermodynamic Properties 
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VAN'T HOFF PLOTS FOR REVERSIBLE FLUIDS 


Appendi y 
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